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We present a simple derivation of the energy formula found by Tan, relative to the single channel 
hamiltonian relevant for ultracold Fermi gases. This derivation is generalized to particles with 
different masses, to arbitrary mixtures, and to two-dimensional space. We show how, in a field 
theoretical approach, the 1/fc 4 tail in the momentum distribution and the energy formula arise in 
a natural way. As a specific example, we consider quantitative calculations of the energy, from 
different formulas within the ladder diagrams approximation in the normal state. The comparison 
of the results provides an indication on the quality of the approximation. 



I. INTRODUCTION 



The remarkable progress in the field of ultracold atomic gases has provided access to a number of systems which 
may display quite new physical properties. One of the most striking example is provided by ultracold fermionic gases 
[l[ and the BEC-BCS crossover. While the physics of Bose-Einstein condensates (BEC) is known from superfluid 4 He 
(and from ultracold bosonic atomic gases), and the one of BCS condensates from superconductors, fermionic gases 
provide systems which display a smooth continuous evolution between these two extremes. This is made possible by 
the existence of Fcshbach resonances which allow to control the value of the scattering length a merely by changing 
the applied magnetic field. For ultracold gases the kinetic energy of the atoms is so small that only s-wave scattering is 
relevant, and it is fully characterized by the scattering length. Since Pauli principle forbids s-wave scattering between 
identical atoms, the scattering length corresponds in most of the experiments performed so far to scattering between 
atoms in different hypcrfine states of a same element. These hyperfine states are often called for convenience "spin 
up" and "spin down". In the case of a wide Feshbach resonance (as it occurs for example for 6 Li and 40 K), the 
closed channel responsible for the Feshbach resonance may be omitted and the system is described by a single open 
channel hamiltonian, where the interaction is characterized by the single parameter a. With the densities produced 
in experiments this interaction has quite a small range compared to the interatomic distance. This single channel 
hamiltonian is particularly interesting because it is at the same time very simple and highly non trivial, and moreover 
it is realized with an excellent precision in these ultracold fermionic gases. Hence we may hope, in getting full control 
and understanding of this hamiltonian, to gain knowledge which may be applicable to more complex hamiltonians, 
such as those encountered in condensed matter physics or in quark matter [2|. 

In this context a simple general expression for the energy of a system described by this hamiltonian has been 
found by Tan Q, which involves only the momentum distribution n a {k) of the particles together with their large 
k behaviour. However the details raise unanswered mathematical difficulties. Nevertheless the expression can be 
checked in limiting cases (see Appendix |A1 for details). This expression is of high interest since it is quite non trivial. 
It is directly related to the well-known problem that, for a contact interaction, the kinetic energy presents a formal 
divergence, because the momentum distribution n a {k) behaves as 1/fc 4 for large momentum (as found for example 
in the perturbative calculation of Belyakov Q), a feature merely linked to two-body physics. Naturally there is no 
divergence in the energy itself, because the interaction energy comes in to compensate this feature from the kinetic 
energy. This is easily seen from a simple finite range interaction model (for example the square well potential), with 
range r$. The momentum distribution decreases more rapidly than 1/k , for momentum beyond 1/Vo which acts as 
a cut-off in momentum space. The apparently singular situation is found when one lets ro — » 0, which should be 
taken as the definition of the contact interaction. In the formula found by Tan, the total energy appears essentially 
as the kinetic energy with this divergence problem removed. Indeed the high momentum part, responsible for the 
divergence, is subtracted out and an additional explicit contribution which contains the scattering length appears. 

This matter has been taken up recently by Braaten and Platter Q who made use of the operator product expansion 
developped by K. G. Wilson. They have in this way derived the energy relation in a compact and formal way. In 
contrast with Tan approach the interaction comes in explicitly. They have also obtained the expression of the known 
@ adiabatic relation for this contact potential hamiltonian. This convenient expression in terms of the scattering 
length and of the coefficient of the 1/fc 4 tail in the momentum distribution has been pointed out by Tan 0] • Quite 
recently it has been used Q to obtain the number of closed-channel molecules in the two-channel model. Also very 
recently it has been derived and studied in details by Zhang and Leggett Q. 

In this paper we present first a derivation of this energy formula, which is simple, explicit, fairly short and avoids 
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unnecessary and uncontrolled mathematical complications. It is in line with Tan approach, in that it deals only with 
the kinetic energy. This is done in a careful way in order to avoid divergences. Our treatment is similar in spirit 
to Rcf. 9]. This simple proof allows straightforward generalizations to more complicated situations, with unequal 
masses, several kinds of particles and two-dimensional space. Next we show how, in a field theoretical approach, the 
1/fc 4 tail in the momentum distribution and the energy formula arise in a natural way. Finally, as a specific example, 
we consider also quantitative calculations within the ladder diagrams approximation in the normal state. There are 
different ways to obtain the energy and we can compare the various approximate results, which gives an indication 
on the quality of the approximation. 



II. DETAILED DERIVATION 



In this paragraph we consider the case where we have only two kinds of particles which are ultracold fermions. The 
volume of the system is assumed to be unity. Namely we have n-f particles with mass m| and ny particles with mass 
TOj. We consider directly the case where the masses rcif and are different since it does not make any problem. 
The positions of the | particles are denoted r^, while those of the J, particles are pj. For these ultracold fermions 
only s-wave scattering has to be taken into account. Hence only interactions between | and J. particles have to be 
considered. Generalizations are considered in the next section. The Hamiltonian reads H = H c + V with (we take 
h=l): 

H c = VA r . VA„, (1) 

2m t ^ • 2mi ^ Pj K ' 

1 i=i 1 j=i 

and 

V = E V ^ - Pi) ( 2 ) 

We assume the interaction potential V(r) to be short-range. For clarity and simplicity we assume that it has a definite 
range r$ and satisfies V(r) = for r > ro- However, just as in standard scattering theory [f|, the results will hold for 
physical short-range potentials where the interaction decreases rapidly enough with interparticlc distance. 

We consider an eigenfunction |$) of this Hamiltonian, having for example in mind the ground state wavefunction. 
However this is not necessary and we can as well consider excited states, which leads to the extension of the results to 
non zero temperature, as pointed out by Tan Q. Let $({ri}, {Pj}) be the corresponding wavefunction, with proper 
symmctrization and normalization to unity. We want to calculate the energy corresponding to this state: 

E = = J dr t d Pj $*({rj, { Pj }) (H c + V) $({r 4 }, { Pj }) (3) 

where dvi dpj is for Y\i j d-Vi dp, with i = 1, • ■ • , w\ and j = 1, ■ • ■ ,n±. 

We follow the basic principle put forward by Tan In Eq. the potential V(r) comes in only in the regions of 
the integration domain where |r m — p n \ < ro, where m = 1, • • • , nj and n = 1, • • • , n^. The overlap between these 
regions is negligible for dilute gases where ro is small compared with the mean interparticle distance d. There are 
n j ni sucrl regions. Compared to the total integration domain, the domain of these regions is of order (ro/e?) 3 . For 
equal populations rif = nj. = kp/6n 2 , this would be of order (fcpro) 3 . Since we have = E<& also inside these 
regions, and that <f> inside these regions has no singular behaviour (this can be checked explicitly since this behaviour 
is obtained from a one-body Schrodinger equation, as we use below), we may neglect the contribution of these regions. 
Actually this has to be done for consistency since @ keeping only s-wave scattering for ultracold gases implies that 
terms of order (fci?ro) 2 , corresponding to higher angular momenta, are neglected. Hence the energy is simply obtained 
by calculating the kinetic energy outside these regions: 

E = f dr t dp, $*({ ri }, { Pj })HM{ ri }, { Pj }) (4) 

J \r m —p n \>r 

We introduce now the Fourier transform i^({k;}, {q^}) of <&({r 2 ;}, {Pj}) with respect to all variables: 

*({r i },{p j })= J2 k '' " >: " " -{M.{qJ; (5) 

{k.},{q 3 } 
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where we use the notation ^ k = (2tt) 3 J dk. The one-particle density distributions n-[(k) and n±(k) are given by: 



n T (fcO = n T Y, M{kJ,{q,})| 2 
{k,^i},{ qj } 

Mn) = H Y lv({ki}:{qj})l 2 

{k»},{q 3 -,j¥l} 



(6) 



This leads to: 



E = y 



(7) 



{k,},{ qj } 



2m T ^ M 2m i ^ 



9iv 



Af=l 



7V=1 



^({ki},{cb-}) 



|r m — p„|>ro 



If we did not have the restrictions \r m — p n \ > tq in the last integral, it would give a factor proportional to Yii <^(ki — 
k'i) Oj ^(^j — an d we would find the standard expression for the kinetic energy: 



Ec = 7T— Y fc i"T( fc i) + 7T— Y 

1 kj 1 qi 



(8) 



which displays a divergence for large wave vector. Hence the restriction is crucial in order to avoid the divergence. 
However, since we want to have the standard expression for the kinetic energy appearing, we will write: 



/ dr l dpj = / dr t dp } ~Y/ dri ^Pj 

J\r m -p n \>r J m n J\r m -p„\<r 



(9) 



where for simplicity we have not written explicitely the exponential integrand appearing in Eq.([7]). As long as all 
the fcj and qj integrations are not performed in Eq.([7]), no singularity appears so we can handle the various terms 
in Eq.@ separately. As we have just mentionned the first one leads to the standard kinetic energy expression, so 
we have just to handle the other terms quite carefully since the limit tq — > is singular. Note that, in the regions 
|r m — p n \ < To, we will deal with the analytic continuation of the expression of the wavefunction for |r m — p n \ > ?*o 
which we will discuss below. This is naturally quite different from its physical value in these regions, since, for these 
last one, the interaction V(r) has to be taken into account. 

Let us consider the term m = n = 1 in Eq.©. In Eq.([7]) only the terms {l/2m{)k\ + (l/2m[)q 2 in the bracket 
will be relevant, and for now on we consider only these ones. Basically we will deal carefully with this two-body 
problem. Accordingly it is convenient to think, as long as we are not back to the many-body problem, that only 
these two particles scatter. Equivalently we could handle all the {m, n} terms simultaneously, but this is much more 
cumbersome presentation. Since they are more convenient to express the condition |ri — p-J < rg, we introduce the 
relative position r = ri — p 1 and the center of mass position R = (m^ri + m^p^/M, together with their conjugate 
variables the total momentum K = ki +qi and the relative momentum k = (mj_kj — m^qx)/M with M = m-j- +mj. 
We have: 



1 



2m 



1 



2m 



9 

■Qi 



^ 2 
2/j 



1 

2M J 



:K 2 



(10) 



where fi = m-[mi/(m-\ + mjj is the reduced mass. The difficulties arise from the k 2 /(2/i) term in the large k limit. 

It is convenient for the discussion to introduce the Fourier transform </?(k, K, {k 2 ...}, {q 2 ...}) of the wavefunction 
with respect to the variables r and R, instead of ri and p x . It is merely given by: 



£(k, K, {k 2 ...}, {q 2 ...}) = ^(k + -Ik, k + -±K, {k 2 ...}, {q 2 ...}) 



M 

Making use, inEq.© andEq.©, of (k 1 -k / 1 ).ri+(q 1 -q' 1 ).p 1 = 
and dki dqi = dkdK., we find that, for the fc 2 /2/i term in Eq. 



(k- 



M 

k').r+(K-K').R , together with dx\ dp x 
we have to deal with: 



Y 5>*(k',K,kq 2 )fc 2 ^(k,K,kq 2 ) f dre l ( k - k ') r 

K,kq 2 k,k' Jr<r ° 



(11) 



drdR 



(12) 
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where kq2 is a short-hand for {k2, ■ • ■ }, {q2, ■ • • }• We have used the fact that integration over all unrestricted position 
variables gives corresponding S functions for the corresponding wavevector variables, as we have mentionned above. 
The problem for large k in this expression is directly linked to the behaviour of the wavefunction for small relative 
distance r, which we consider now. 

When the relative distance r is small compared to the mean interparticle distance d, the dependence of the wave- 
function <J?({i"i}, {Pj}) on r is given by the solution of the relative motion of the two-body problem. For ultracold gases 
the energy corresponding to this motion is nearly zero and the corresponding wavefunction is, for r > ro, proportional 
to ip{ r ) with ip(r) = l/ r — l/ a where a is the scattering length This form is actually valid provided r is small 
compared to a typical particle wavelength, that is r <C d. Since we have ro <§C d, there is a range of validity for this 
form for r > r . We are also interested in the case of large positive a, where a two-body bound state exists with 
wavefunction proportional to exp(— r/a)/r. In this case the above form requires r <C a. There is again a range of 
validity for this condition for r > r$ since we are in practice interested in physical situations where the scattering 
length a is large compared to the potential range ro. Hence for small r we have: 

$({r,}, { Pi }) = lft(r)*(R, {r 2 ...}, {p 2 ...}) (13) 

Since in the evaluation of the J r<r term of Eq.©, we will let r — -> 0, only large values of k are relevant. Indeed if 
we considered only bounded values k < k c , this term would go to zero for ro — » (and the k integral would converge 
because of the cut-off k c ). Hence we have to consider k — > oo. In this case the dependence of (p(k, K, {k 2 ...}, {q 2 ...}) 
on k is entirely linked to the short distance behaviour on r of $({ri}, {Pj}) given by Eq. lfTB")) . We have in this limit: 

$j(k, K, {k 2 ...}, {q 2 ...}) = VfWMK, {k 2 ...}, {q 2 ...}) (14) 

where tpF(k) and $_f(K, {k 2 ...}, {q 2 ...}) are the Fourier transform of i/j(r) and $(R, {r 2 ...}, {p 2 ...}) respectively. 

As it is well known, the Fourier transform of 1/r is An/k 2 . On the other hand, for the calculation of 
k 2 y>(k, K, {k 2 ...}, {q 2 ...}) in expression (fTSj) . we do not have to take into account the constant —1/a in the wave- 
function since it gives zero when we apply on it the kinetic energy operator A r (corresponding to the factor k 2 ). All 
this Fourier transform calculation can be done quite carefully by multiplying the wavefunction ip(r) by a convergence 
factor e~ vr , then letting 77 — + . This confirms the above results. Finally, when handling expression (|12p. we have 
merely k 2 p(k, K, {k 2 ...}, {q 2 ...}) = 4^$ F (K, {k 2 ...}, {q 2 ...}). 

Physically: 

p(k)= l^(k,K,kq 2 )| 2 (15) 

K,kq 2 

is the isotropic probability distribution of wavevector k in the relative motion of the two f and J, particles we are 
considering. We have found that its leading behaviour for k — > 00 is p4/k 4 with: 

Pi = (An) 2 J2 l*F(K,kq 2 )| 2 (16) 

K,kq 2 

On the other hand, with respect to k' integration in expression (I12p . it is more convenient to go back to r space 
through e~ lk r ^'J(fc') = i/j(r), which gives a factor tp(r)<&* F (K, {k 2 ...}, {q 2 ...}). Since ip(r) is isotropic, we can 
perform cxplicitely the angular r integration J df2 r e tkr = Air sin(fcr) / (kr) . Finally, for expression (|12p . we are left 
with the calculation of dr r sin(kr)ip(r). This r integration is easily performed, leading to: 

k[ drrsin(fcr) ( - - - ) = 1 - ( (1 - — )cos(fcr ) + ^-sm(kr ) ) = 1 - f(k,r ) (17) 
\r a J \ a ka 1 



We write now, with this result, the partial contribution Enk to the energy Eq.([7]), coming from the term m = n = 1 
in Eq.@ and where we retain only the k 2 /2fj, term in Eg. ([TO)) . We obtain: 

E^ = --J^^ (i_ /(fc>n) )) = -_y o dk P4 + £L-l dkf { k,r ) (18) 

In the first term in the right-hand side, we have put a cut-off k c since this integral diverges when k c — ► 00. This 
divergence will just compensate in the final result the above mentionned divergence in Eq.©. Hence the global result 
will be convergent as expected. 

If in the last integral of Eq. (fT8"|) we were setting r = 0, we would get f(k,r ) = 1 and a divergent result. Instead 
this integral is perfectly convergent when we calculate it explicitly for ro 7^ and then take properly the ro — > limit. 
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Indeed we see that it is not divergent because f(k, ro) involves oscillatory functions. The cos(fcro) term merely gives a 
result proportional to <5(r ), where S(x) is the Dirac distribution. Since r ^ the contribution of this cos(fcr ) term 
is zero. On the other hand, owing to: 



dx 



sin 2! 



(19) 



the sin(fcro) term gives a contribution 7r/2a to the integral, leading to a contribution p4/(87r/ia) to En^. As we 
have mcntionncd above, we could improve the presentation of our handling of the k integration, to deal with perfectly 
defined integrals, by introducing a convergence factor exp(— rjk) and then let r/ — > 0+. Physically this would correspond 
to regularize the wavefunction in the r — > limit. This would confirm our above results. 

We have not yet taken into account the K 2 /2M term in Eg . (flU)) . because there is no singular behaviour associated 
with it. The corresponding contribution E\\k of the term m = n = 1 in Eq.([9|) goes to zero as ro — » 0, as it is obvious 
directly and can be checked by following the same procedure as above. Finally one sees easily that, for all the other 
kinetic energy terms (l/2mj) J2m=2 + 0-/^ m l) T^=2 In m the bracket of Eq. (JJJ) , the contribution from the term 
m = n = 1 in Eq.® also goes to zero when ro — > 0. Hence the total contribution En of the m = n = 1 term is 
merely E u = E llk . 

We rewrite now the sum En in terms of the variables corresponding to particle | and j, instead of the relative and 
center of mass variables. We have: 



-En — — — 



2u ^ 

^ k,K 



f 4 (K) 

k 2 



P4 

8nfia 



with again a cut-off k c understood for the summation over k, and where we have introduced: 



P 4 (K) = (47r) 2 ^|$ F (K,kq 2 )| 2 

kq 2 

related to p 4 by p 4 = P 4 (K) . With 1/7* = \jm\ + 1/m;, we have: 

E 1 V ^(K) 1 V P 4 (K) 



k,K 



k,K 



k 2 



Pa 
87r/ia 



(20) 



(21) 



(22) 



We make use of Eg. tfTTJ)) to go back to the ki and qi variables. In the first term we change the summation variables 
from {K, k} to {K, ki = k + (m^/M)K}, making use of dkdK = dk± dK. In particular we have 1/fc 2 = l/(ki — 
(m-f/M)K) 2 . However we can use the identity: 



1 



1 



(r-A)* 







where A is any fixed vector (this result is obvious by changing r into — r), to write, X^kilV^ 2 — V^i] 
K. Proceeding in the same way with the second term of Eg. (f2"2")l . we obtain: 



F = 1 yft 1 yft 1 

^ n 2m* ^ i- 2 O.m , „? T 



1"* 



P4 



(23) 
at fixed 

(24) 



with a cut-off k c understood for the summation over ki and qi. 

In conclusion we have calculated all the contributions coming from the presence of the term m = n = 1 in Eq.®. 
These are just the three terms, proportional to p 4 , in Eq. (|24l) . To summarize, for the expression of the energy, we are 
back to Eq.0 with the restriction |ri — p\\ > ro removed, and the above three p 4 terms added. 

We have just to repeat the same argument for all the other restrictions |r< — Pj\ > r$. In this way we obtain rifrij 
analogous p± terms, which are naturally all equivalent after a change of variables. Taking Eq.® into account, this 
leads to the final expression for the energy: 



E 



1 



2mt ^-^ 

1 L- 



\k 2 n^{k) 



"4 
k 2 



+ 



2m 



E 



Q 2 ni(q) 



«4 
9" 

9 



7l 4 

8irfia 



(25) 



where n 4 = ri|njp 4 = lim^oo fc 4 ri|(fc) = lim^oo k A n^(k). Indeed, from Eq.®, rtf(fc) and ni(k) behave in this way 
for large k since for example ki scatters with all the q^, and each scattering brings a contribution p^/k 4 . That is any 
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of the rij particles scatters with any of the particles. Hence the summations in Eq. (|25|) arc perfectly convergent. 
This formula is the simple generalization of the formula found by Tan Q to the case where the two species of involved 
particles have different masses. 

In this derivation we introduced a cut-off to manipulate separately each contribution. This makes an easier presenta- 
tion for the derivation. However this is just a convenience. We could avoid it by handling all the terms simultaneously. 
The presentation would be much awkward, but we would only deal with well defined convergent integrals, without 
any need for a cut-off. 



III. SIMPLE GENERALIZATIONS 



It is first worthwhile to note that the above derivation did not make use of the statistics of the particles. Hence 
the result is valid for bosons as well as for fermions. Naturally it is artificial for bosons to consider only scattering 
between different species, although it might just happen that this scattering is the dominant one. Nevertheless for 
bosons there is no reason to exclude scattering between particles belonging to the same species. If we consider first 
the case of a single bosonic species, with n particles of mass m in the unit volume and with scattering length a, we 
can follow the same procedure as in the preceding section. We can write the equivalent of Eq.([7]) for a single species, 
with a restriction |r, — Vj\ > r a working now between any of the n(n — l)/2 couples of particles. Taking care of the 
restrictions in the same way as in the preceding section, we end up with: 

E = ^-T.\ k2n ( k )- ! n}+^- ( 26 ) 

2m ^ L k 2 \ 87rma 

k 

with n 4 = n(n — l)p 4 = lim^^ k n(k). 

We can then generalize this result to any mixture of N boson species with m particles of mass m, in the unit volume 
(i = I , • • • , N) and interaction between species i and j characterized by scattering lengths ay. We find in the same 
way: 

N 1 1 

E = Y— YWn^k)-^] +— Y (27) 
^ 2rrii ^ L W k 2 J 16tt j-< V ; 

where fiij = miirij / (mi + rrij) is the reduced mass for the ij scattering, na = linifc^oo k A rii{k) and na = riuj- 
Naturally the formula is more complex for these mixtures, since it requires the knowledge of the N(N + 1)/2 constants 
n 4ij = n 4ji associated with the i — j scattering. The same result works for ultracold fermionic mixtures, except that 
we have to set n^u = because of Pauli principle. 



IV. THE 2D CASE 



While in ID situations the kinetic energy converges, and indeed in exact solutions of many-body problems the 
energy is precisely calculated by evaluating the kinetic energy, there is in 2D a divergence analogous to the 3D case, 
except that the divergence is logarithmic in the 2D case. We show here how the procedure followed in 3D can be 
extended to this 2D case. Actually there is not so much difference since the space dimensionality does not appear 
in the principle of the procedure. The changes appear only when one comes to pratical matters. We keep the same 
notations for the variables, but naturally we have to deal now with two dimensional variables and integrations. 

First the expression of the wavefunction ip(r) for the relative motion at small distance r > r$ is modified into 
ip(r) = ln(a/r), which is solution of the 2D equation Aip(r) = 0. The length a for which ip(a) = plays the role 
of the scattering length and is naturally obtained from the interaction potential. The Fourier transform of ip( r ) IS 
ipF(k) = 27r/fc 2 , omitting the irrelevant Fourier transform of In a. Similarly the large k behaviour of p(k) is p4/k 4 , 
with (27r) 2 appearing in Eg. (fTB)) instead of (47r) 2 . Then the 2D angular integration gives f Q * e lkr = 2ir Jo(fcr) where 
Jo(x) is the first kind Bessel function. Hence, instead of Eq. lfTS)) . we obtain: 

Pi f dk f r ° a , N p 4 [ K dk f krn ka , , 

E ^—yJ Jo drrXn -r Ukr) = -^kL tI dxxln —Mx) (28) 

where in the last expression we have changed to the variable x = fcr, and we have naturally to take the cut-off k c — > oo 
at the end of the calculation, when the standard kinetic energy term is included, as we have done in section ITT] The 
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logarithmic divergence arising in this standard term is expected to be compensated by the contributions of all the 
terms similar to Eq. (|28p . Indeed, integrating by parts, we have: 



dk 



kra 



dxx In a; Jo (a;) 



Ink 



kr 



dxx In a; Jo(x) 



dkk In k ln(fcro) Jo(fcro) 



(29) 



To evaluate the first term we use |LT| (when the integral is not absolutely convergent, we define it as above with an 
exponential convergence factor, with an extremely weak decreasing behaviour): 



f 

Jo 



dx x^ Jq (x) 



r(i^) 



(30) 



(where T(x) is the standard Eulcr gamma function) together with \nx = lim e ^o(x e — l)/e. Note in particular that 
/ °°dxx Jo(x) = 0. This gives J °° dxx In x Jo (x) = —1, which is also consistent with the fact that the Fourier 
transform of lnr is — {2ir)/k 2 . Hence we find indeed that Eq. (|29p provides the required compensating term to avoid a 
divergent result. The other terms lead to constants, so we may take immediately in their expression the limit k c — > oo. 
First, except for the prefactor — p4 / (Air n) , the factor of In a in the second term of the integral in Eq. (|2"5)l is: 



dk 

T 



kr 



dx X Jq (x) 



dy 

o v Jo 



dxx Jq(x) 



dyy lny J (y) = 1 



(31) 



by the change of variable y = kr$ and by integrating by parts. This means that the divergent contribution of the 
second term in Eq. (|28[) is proportional to ln(fc c a), as could be expected from dimensional analysis. 
Finally in Eq. (|28|) we have still a contribution proportional to In k : 



f°° dk f kr ° r 2 f°° 

I —— In k dx x Jo (x) = — dk k In k Jo 
Jo k Jo 2 Jq 



(kr ) 



(32) 



again by integrating by parts. Gathering this term and the second term in Eq. (|29jl . we have: 

„2 



dk k \nk\n(kro) Jo(kro) + / dk k In 2 k Jo(fcro) 



d(kr ) (kr ) In (kr ) J (kr ) (33) 



where, in the last step for the term proportional to ln 2 ro, we have used again J °°dxxJo(x) = 0. The remaining 
integral is calculated again by making use of Eq.lUP and is found to be (1/2) J °°dxx In 2 x J (x) = C-hx2 ~ -0.116, 
where C = 0.577216 • • • is the Euler constant. Gathering the above results we obtain: 



P4 w , v , In2-C 
Euk = —. — m(fc c a) + P4 — 

47T/i 47T/i 



(34) 



Finally, just as for the 3D case, we want to go back to the variables corresponding to particle f and |. It may be 
seen that this does not lead to additional contributions. This leads us finally to: 



E 



1 



2m-f fc c 



lim 



'fc<fc c 

^ k 2 ni(k) - — ln(fc c a) 

. k 



lim 

2l7ll qa^oo 



q<qc 



m ^jr (35) 



where we recall that 714 = n|n^4 = Mm^-toc k^n-\{k) = lim^oo k^n^k). This result can be checked explicitly in the 
molecular case (see Appendix A). 



V. FIELD THEORETIC APPROACH 



Let us see now how the above expression for the energy arises in the field theoretic formalism. To be specific we 
will restrict ourselves for simplicity to the case of major interest, namely the one of two fermionic species with equal 
populations, so that n-\(k) = ni(k) = n(k). Since all the f and j quantities are equal, we do not write explicitely this 
index. We consider also the non zero temperature case since it does not make any problem. We deal first with the 
case of a normal system, and then extend the results to the superfluid case. 
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A. Normal state 



We consider first the large fc dependence of n(k). For free fermions at temperature T, the density distribution has 
an exponential tail proportional to e~ k l 2mT '. However, as pointed out in Ref. interactions modify this behaviour 
and give rise on general grounds to a 1/fc 4 dependence which dominates the exponential tail. This is explicit in the 
weak coupling domain where the interaction can be treated perturbatively and the scattering length a is small, as it 
has been done by Belyakov [H at zero temperature: 

/ 2 \ 2 k 4 

n(k) = (^-k F a) ^ (36) 

where kp is the Fermi momentum = ?i| = k F /(6ir 2 ). 

In the general case the distribution n(k) is obtained from the temperature Green's function G(k, icu n ), where 
u> n = (2n + 1)ttT (with n an integer and kg — 1) is the Matsubara frequency, by: 

n(k) = T ^ G(k, iu n ) e^" T (37) 

n 

where r — ► 0+. We separate out in this equation the free particle contribution by writing Dysons 's equation: 

G(k, iw n ) = G (k, iuj n ) + G (k, io; n )E(k, iw„)G(k, iw n ) (38) 

where Go is the free particle Green's function Go(k, iuj n ) = [iu> n — + fi] , with = k 2 /2m the free particle kinetic 
energy, and E(k, ioj n ) is the self-energy. The first term in Eq. (j38|) gives in Ea.([37)l the free particle contribution, 
namely the Fermi distribution. For large k its exponential tail mcntionncd above is completely dominated by the 
algebraic decay 1/fc 4 that we will obtain. Hence we are left only with the second term. For large k, implying a 
large kinetic energy for the particle, we expect the effect of interaction to be small in the same spirit as the Born 
approximation in this regime. Since S(k, iu> n ) describes this effect we expect it to be small. Hence we may to lowest 
order replace in this second term G by Go- On the other hand we replace, in a standard way fl"ot | . the summation over 
Matsubara frequencies by a frequency integration over a contour C encircling the imaginary axis in the anticlockwise 
direction. This leads to the following expression for the dominant contribution to n(k) at large k: 

«k) = -± [*,/£) -. £(k '"> (39) 
2«r Jc T (cj-e k +fj,) z 

where f(x) = l/(e x + 1) is the Fermi distribution function. The contour can be deformed into the sum of a contour 
enclosing the positive frequency Rew > half-plane and another contour enclosing the negative frequency Rew < 
half-plane, both being in the clockwise direction. Closing the contour at infinity in the Rew > half-plane is allowed 
by the presence of the Fermi distribution f(u)/T), and closing it at infinity in the Re a; < half-plane by the presence 
of the e WT factor, which is then omitted since it does not play any other role. 

Now contributions from the double pole at w = f t - ^ will contain from the Fermi distribution a factor e~ k l 2mT 
which makes them negligible. Similarly E(k, lS) has also singularities with frequencies which arc large and positive, 
when k is large. The Fermi distribution f(uj/T) will again make their contribution exponentially small. On the other 
hand, as it is shown in details in Appendix [Bl E(k, u>) has a pole at u) ~ — e k for large k, that is deep in the Re u> < 
half-plane. Roughly speaking this pole appears since a t particle with large k and lu will scatter with a J, particle 
with essentially opposite parameters — k and —u>. This happens because, just as in section HH the large values occur 
for the relative motion, but not for the center of mass motion. Then the J, particle propagator has a pole at u ~ — e^, 
which produces a pole for E(k, uS) at the same frequency. 

For this pole, we have merely (lu — e k + /i) 2 — (2efe) 2 = k 4 /m 2 and /(— e&/T) ~ 1. If we call Jfe the corresponding 
residue of E(k, uS), we obtain in this way: 

n(k) ~ m 2 R s -1 (40) 

Hence we see that the 1/fc 4 dependence of n(k) at large k emerges quite naturally in this approach. We note that 
our analysis may seem inconsistent, since we have first argued that E(k, u>) should be small for large fc, but then 
considered a pole of E(k, u>), in the vicinity of which it is large. The justification is that, as long as the contour is far 
away from the pole, E(k, u) is indeed small everywhere on the contour. Making use of the pole is then a convenient 
way to calculate the contour integral. 
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We consider now the expression of the energy. It can be written in terms of the Green's function by a standard 
formula [To| which bears from the start a strong analogy with the expression found by Tan Q : 



E = T^ fan + M + £ fc) G(k, iuj n ) j 



(41) 



where again r — > 0+. This equation includes an overall factor 2, coming from summation over spin. The essence of 
this formula is that, calculating ioj n G(k, iu> n ) by Hciscnberg equations of motion, one finds E c + 2Ei nt , that is the 
kinetic energy E c plus twice the interaction energy E lni . Adding the kinetic energy E c , which gives the Ck term in 
Eq. (|4"Tj) . leads to 2E = 2(E C + Ei nt ). In practice Eq. (|4"Tj) is not so useful since one has to deal carefully with the 
divergent behaviour occuring for large uj n and large k, which is quite painful numerically. In this respect the number 
equation Eq. (|37p is much more convenient. 

Now, to single out clearly the kinetic energy contribution, we can rewrite Eq. (|4ip as: 

(42) 



E = 2TJ2 G(k, iu n ) e l ^ T + TJ2 («"« + M - e fc ) G{k,ito n ) e 1 ' 



The first term is just the kinetic energy E c = 2 ^ k efc rt(/c). In the second term, which is just -Eint, we can again 
replace the Matsubara frequency summation by the same contour integrals as above. Moreover the quantity appearing 
in this summation is just Gq uj)G(k, lu) = 1 + S(k, u))G(\l, uf). However the term 1 gives a zero contribution to 
these closed contour integrals. Hence we are left with: 



E = E r 



—y 

k 



dw/(y)E(k,w)G(k J w) 



(43) 



We can now analyze the behaviour of the interaction energy term for large k, as we have done above for the particle 
number. In this range E(k, u>) is small so we may replace G(k,u>) by G?o(k, ui). Hence we have to deal with the same 
expression as in Eq. (|39p . except that Go(k,uj) is not squared. Following the same arguments we obtain for large k: 



/ d W /(^)E(k,w)G(k,a;) 
2m J c T 



(44) 



Hence we have shown explicitely that the large k behaviour of the interaction energy comes in to cancel the divergent 
behaviour of the kinetic energy. Naturally this is the expected result. Accordingly we may write: 



E 



2m 



n(k) 



2k 2 



k 



-L / dw/(^)E(k,w)G(k > a;) 

lVK J c T 



(45) 



where, from Eq. (|4"0|) and Eq. (|4"4"l) , both brackets give convergent integrals for large k. The first term is just the kinetic 
energy from which, according to Eq. (|4"D|) . the large k behaviour has been subtracted. By comparison with Eq. (f2"5)l . 
the second term in Eq. (|45[) is just mR-^/ (Aita), a result by no means obvious. 



B. Numerical calculations in the ladder approximation 



We are now in a situation where we can calculate the energy from three different formulas, as we explain just below. 
For the exact theory they would naturally give the same result. However if we use an approximate scheme, the results 
are expected to be different. The differences can, in some rough way, be seen as a measure of the errors resulting from 
the approximation. Hence it is of interest to make such a comparison. Specifically we have taken as an example the 
ladder approximation (l3l . [l4| at unitarity and we have performed numerical calculations of the energy as a function 
of temperature in this approximation. 

We use three different ways to obtain the energy. Our first expression for the energy is the general one Eqs. (|43l45p . 
The second one is the energy formula found by Tan, that is Eq. (|25|) with equal masses. Finally we use the non- 
interacting gas scaling relation E = — (3/2) fi between the grand potential f2 and the energy E, which is also valid at 
unitarity [15| due to the lack of any microscopic energy scale. We calculate the grand potential itself via the density 
?i(/i, T) (obtained by integration of n(k)) using the thermodynamic relation f2(/x, T)/V = — dfx' n(// , T) (we have 
used the fact that f2 goes to in the high temperature regime fi/T — > — oo). As noted above this last method is 
much more convenient numerically than the two others, since there are no convergence problems for large k in the k 
summation. 
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The results of these numerical calculations are shown in FigfT] We have checked that we recover the proper high 
temperature behavior (virial expansion) with the three methods, as it should be. This common evolution toward 
the virial result is clearer in the inset, where the results are shown with an extended scale, but this convergence is 
fairly slow. The results of the three methods are fairly close, which is coherent with the fact that we expect the 
ladder approximation to be reasonably good, mostly toward higher temperatures. Indeed at lower temperature the 
differences between the three results are growing, mostly for the result of the Tan formula. This is clearly due to the 
proximity of the superfluid critical temperature, which occurs at T/Tp — 0.243 in this approximation. 

The major numerical problem with the use of Eq. (|25p is to obtain the proper coefficient of 1/fc 4 in the large 
momentum behaviour of n{k). This requires in particular to be sure that the asymptotic regime is reached numerically. 
The same problem arises in Eq.([45p since one must find numerically that, for large k, the kinetic energy is exactly 
balanced by the interaction energy. In the present case of the ladder approximation, we can check that the asymptotic 
regime as been reached since we have analytical expressions. In order to check numerically this large momentum 
behaviour of the contributions coming into the two brackets of Eq. (|45p . we have performed, as we indicate specifically 
below, an expansion of the occupation number n(k) and of the interaction energy term up to order k~ 6 and fc~ 4 
respectively. 

As indicated in the above subsection the self-energy is small in the large momentum and frequency limit and we 
can expand the Green's function in powers of the self-energy. Hence we have for the occupation number and the 
interaction energy term: 

n ( k )=-^-J c duj f{^) (Go(k^) + Go(k, w )S(k, w )G (k,^)+0(S] 2 )) (46) 

~inj c dujf (I) E (k ' U) ° (k ' U) = ~ki Jc (?) (S (k ' U) G ° (k ' Uj)+ ° (E2)) (47) 

We start from equation (Bl) and following the arguments given above and detailed in Appendix B, we neglect the 
contribution coming from the pole of the Green's function, that gives exponentially small terms, and we only consider 
the contribution from the cut of the vertex Y (K, 0) on the real frequency axis. Following [l4j|, we call this self-energy 
contribution Sr: 

S r (k, u) = - J dnb (^pj Imr(K,Q + ie) G (K - k, - w) (48) 

with Q min = |^ - 2(i, and e -> 0+. 

The Bose distribution effectively limits the frequency f2 to be at most a few T and the JC-integral is therefore 
effectively bounded. Similarly to what is done in Appendix B, we consider the singularities of the self-energy located 
at u> f=a — £k- In the large momentum k and large frequency u> limit, we expand the free Green's function in Eq. (|48p 
in power of (oj + £k + fi)~ , leading to: 



E r ( fc ,a,) = V ° A l\ p+1 (49) 



where 



(50) 



C v (k)= V / dnb I - ImT (K, fi + ie) x + il + 2/i 

7T ^ Jn min \TJ \ m 2m 

In this limit, the contribution to the self-energy due to the cut of the vertex T (K, £1) is a sum of poles of order p+ 1 
located at u> = — £k — /i. In Eq. (|46[) . only the second term is of interest, since the first one gives the Fermi distribution 
which does not contribute to the algebraic tail. So we use the expansion Eq. ([49|) in Eq. ([46|) and find the following 
expansion for large wave vector k: 




11 



0.8 



0.6 



£ 0.4 



0.2 





0.2 0.4 0.6 

T/T c 

F 

FIG. 1: Energy of the gas at unitarity as a function of reduced temperature, obtained from Eq. (|45[l (dashed line), Eq. (|25p 
(dashed-dotted line) and from the scaling relation (full line) E — — (3/2) Q, with the grand potential Q obtained by integrating 
the particle density. The inset shows the results with an extended scale. 

The first three terms of this expansion give the dominant contribution. Indeed, it is easily checked that the term 
of order E 2 in Eq. ((46|) give contributions at least of order fc~ 8 . Co does not depend on k while the k dependence of 

,2 

C\ vanishes after angular integration. For the coefficient C2 we can write C2 (k) = ^27^ + B 2 . Finally we find for the 
particle distribution in the large momentum limit: 

n (fc) = C -^ + (2d + 3A 2 ) — + f -^J (52) 

Comparing with Eq. (|4"0|) we see that Co is identical to the coefficient i?£ defined in the above subsection. We have 
checked the numerical convergence in the large k limit by comparing the second term in our expansion Eq. (|52p with 
our numerical results. The coefficients are identical with a precision of less than one percent. 

We can perform the same analysis for the interaction energy term in Eq. (|47[) and we find the following expansion : 

" 55 L *" (?) S "> ° * " " C °F " < C ' + A * TF + ° (1?) (53) 

Here also we have checked the numerical convergence with a precision of less than one percent. 

C. Superfluid state 

Compared to the normal state, we have now to modify the Dyson's equation Eq. ([38|) to take into account [TcJ the 
anomalous self-energy A(k, iu) n ) and the anomalous propagator F + (k., iu n ). We have: 

G(k,ioj n ) = G (k,iu! n ) + G (k, iw n )E(k, iw„)G(k, iuu n ) +G (k, iw„)A(k, iu> n )F + (k, iu> n ) (54) 
F+(k,jw n ) = G (-k,-iw n )S(-k,-ia; n )F + (p)-Go(-k,-iw n )A*(k,ia;„)C(k 1 ia; n ) (55) 

Again for large k, we have to lowest order G(k, iuj n ) ~ Go(k, icu n ) and S(k, iu>„) small. For such a large k we expect 
this normal self-energy S(k, iu) n ) to behave as in the normal state, discussed in the preceding subsection. On the 
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other hand the anomalous self-energy is expected to go to a constant, without any singular behaviour. This implies 
from Eq. ([55)l that F + (k, iuj n ) is small. In this case the first term in the right-hand side of Eg. ([55)1 is negligible and 
from the second term we have: 



F + (k 1 iu! n ) ~ -G (-k, -i<j n )A*(k,iuj n )Go(k,iuj n ) 

Substitution into Eq. (|54]) gives: 

G(k,iw„) ~ G Qs.,i(j n ) + Gl(k,iu> n ) [S(k, iu n ) - G (-k, -iu;„)|A(k, iw„)| 2 ] 



(56) 



(57) 



Hence, in addition to the pole at to — e^. coming as discussed above in E(k, itu n ), the last term gives also a pole at 
the same place from the explicit factor Gq(— k, —iu> n ). This leads to the large k behaviour: 



1 



n{k) ~ m 2 [i? s + IAo.1 2 ] 



(58) 



where Aoo = lim/ s _ +00 A(k, — £&). Actually the additional term is well-known in the case of BCS theory which gives 
indeed, with standard notations, Ufc — v 2 — (1/2)(1 — Cfc/^fc) — m 2 A 2 /fc 4 for large k. 

Going now to the expression for the energy, there is naturally no modification in Eq. (|42p . On the other hand we 
have now, in the second term of Eq. (|4"2")) involving Gq 1 G, to make use of the appropriate Dyson's equation Eq. (|Sl 
This leads, instead of Eq. (|15|) . to: 



E = E c -^-Y j f dwf{%) [E(k,c;)G(k, W ) + A(k,^)F + (k,^)] 



(59) 



Now the large k analysis goes essentially as in the normal state, and is similar to the one for the particle distribution 
since there is only a factor Go which is different in the two calculations. This factor leads to a factor — k 2 /m in the 
present case, compared to the particle distribution, that is we have for large k: 



and the generalization of Eq. (|45|) is: 
k 

+ E 



[Ri 



|A n 



k 2 



k 2 

2m 



2fc 2 



(60) 



(61) 



[Re 



k 2 



2^ / dwf(^) [S(k,w)G(k,w)+A(k,Kj)J^(k,iw)] 



VI. CONCLUSION 



In this paper we have investigated the large momentum algebraic tail in the particle distribution and the energy 
formula found by Tan associated with this tail. We have provided a simple derivation of this energy formula, which 
rests on the fact that, in evaluating this energy for the short range potential under consideration, the interaction 
energy contribution is zero in most of phase space, which makes it negligible. Hence only the kinetic energy has to 
be calculated. The basis of the derivation is the careful subtraction of the kinetic energy divergent contribution for 
intcrparticle distance less than the potential range. This derivation is easily generalized to particles with different 
masses, to arbitrary mixtures, and to two-dimensional space. We have then shown how the algebraic tail arises 
naturally in the field theoretical many-body approach, from the analytic structure of the self-energy. Making use 
of these ingredients we have shown how, starting from a standard general expression of the energy in terms of the 
Green's function, one obtains a formula with a kinetic energy part which has the same structure as in the formula 
found by Tan. This has been done both in the normal and in the superfluid state. Finally we have taken the various 
exact formulas allowing to obtain the energy, and we have compared in the normal state at unitarity the resulting 
numerical values obtained within the ladder approximation. 
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APPENDIX A: LIMITING CASES 

In the BEC limit a — > 0+, one is led to consider the problem of a single molecule. The normalized wave function 
is -0(r) = (2na)~ 1 ^ 2 e~ r l a jr and the density distribution is given by the square of its Fourier transform n-t(k) = 



ni(k) = 8wa~ 1 /(k 2 + a~ 2 ) 2 . The coefficient of the fc~ 4 tail is n± = 87ra _1 . One checks that indeed the energy 
E = 2 ^ fe (fc 2 /2m)(n(fc) — 714/A; 4 ) + n^/ (Anma) = — I /(ma 2 ) gives the proper binding energy Et, = l/(mo 2 ). 

In the weak coupling limit a — > 0_ the energy (per unit volume) can be expanded in powers of kpa (the contributions 
which would come from BCS pairing are exponentially negligible). To second order [1 01 ] : 



k 2 

E = n— — 

m 



3 kpa , ., ,n 
ro + ^ + b 2 (k F a) 2 



where n = k F /3ir 2 and 62 = 2(11 — 21og2)/357r 2 ~ 0.0556613. Second order perturbation gives also: 



n T (fc) = ni(k) = 9(k F - k) + (k F a) 2 n {2) (k/k F ) 



(Al) 



(A2) 



where n^(k) is a lengthy analytical expression given in Refs.[J, [T3]. For large k, (k F a) 2 h {2) (k/k F ) ~ n^/k 4 with 
(2k F a/3ir) 2 . Since there is no first order correction to the density distribution, the first order correction to the 



,(2) 



energy, in the energy relation Eq. (|25|) . is merely given by the last explicit term n 4 /4irma = nan 2 /m, which coincides 
indeed with the (mean field) first order term in Eq. (|Al[) . 



To check this relation to second order, we need the coefficient 



.(3) 



of the tail in the third order contribution 



(k F a) 3 h^ 3 \k/k F ) to the density distribution, which is not available. However we can extract it from the adiabatic 
relation 

[EBU "4 = 4nma 2 (dE/da), which gives n\ = 8Trb2nkj 



2m 



(k F a) 2 n (2, (k/k F ) 



which implies: 



dx 



2k F a 
3tt 



(X) 



1 



4 

9^2 



Hence we have to check that: 



2b 2 nk F a 2 k 2 F 2 

£ — = bin— —(k F a) 

m m 



(A3) 



(A4) 



We have checked numerically this equation with a seven digits precision corresponding to all the digits given above 
for b 2 ■ 

Finally in the 2D case, in the BEC limit, one has again a single molecule problem. The normalized wave function 
of the bound state with energy E = — l/(2/ia 2 ) is ip(r) = (a^/ir)" 1 Ko(r/a). For small x the Besscl function 
K (x) ~ \n(2e~ c /x), so that a = 2ae~ c . The density distribution is n^(k) = n±(k) = 47ra 2 /(l + (fca) 2 ) 2 , so that 
714 = 47r/a 2 . With these ingredients Eg. ([55]) is easily checked. 



APPENDIX B: ANALYTIC STRUCTURE OF THE SELF-ENERGY 



We discuss here in details the analytic structure of S(k, uS) in order to justify the assumptions we have made above 
in section IV Al It is useful to consider first the case of the ladder approximation [l(| [H| > where the structure is 
explicit. In this case the self-energy is given by: 

S(k,iw n )=T^r(K,i^)Go(K-k,^-w n ) = ^- I dfifc(p) Er(K,O)G (K-k,f2-io; n ) (Bl) 
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where in the last step we have replaced the summation over bosonic Matsubara frequencies uj v = lixvT ', with v being 
an integer, by an integral over the contour C introduced above, with b(x) = l/(e x — 1) being the Bose distribution 
function. The vertex T(K, CI) is given by: 



r-'<K.«>=£+s: 



T]TGo T (k', lWm )Goi(K - k',0 - iuj m ) - ^~ 



(B2) 



Deforming contour C as we have done above in section IY A[ we will from Eq. (|Bl[) express S(k, uj) in terms of the 
singularities of T(K, CI) and of Go(k, uj), which are on the real frequency axis. On one hand Go(k, uj) has a simple pole 
at u> + fx = ej;. On the other hand the singularities of T(K, CI) correspond first to the continuum of scattering states 
with energy Cl + 2/.i = ey* + eK-k' > K 2 /4m, arising from the product of the two Go in Eq. (|B2[) . In addition there is 
the possibility of a bound state of the two particles, corresponding to a zero of the right-hand side of Eq. (|B2|) . As a 
result the spectrum of these singularities has some lower bound f2 m in- On the other hand for large K their frequencies 
are bounded from below by CI ~ K 2 /4m, corresponding to the kinetic energy of the mass center. 

Now let us first consider the contribution from the pole of Go(K — k, CI — iuj n ) in Eq. (|Bl[) . located at CI = iuj n + 
eK-k — M- The Bose factor b(Cl/T) will produce a factor /(eK-k/T"), which for large fc implies an exponentially small 
factor e~ fe / 2mT . This makes the corresponding contribution to the self-energy irrelevant for the algebraic tail of n^. 
The only way to avoid this is to have also a large K w k, so eje-k is not large. However in this case we will obtain a 
factor r(K, iuj n + £K-k — yu) . The singularities of S(k, uj) we are looking for are obtained from this factor by continuing 
the imaginary frequency iuj n to the real uj axis. However we know that, when K is large, the singularities of T(K, CI) 
are located at large frequencies of order of K 2 /Am. This implies that the corresponding singularities of E(k, w) will 
be for oj + e-K-u — M ~ K 2 /Am ~ k 2 /4m, that is uj ss k 2 /4m since eic-k is not large. As explained below Eq. (f5T]) these 
large uj singularities in S(k, uj) give only exponentially small contributions to and again do not come in the fc -4 
tail we are looking for. Hence we conclude that the pole of Go in Eq. (|Bl[) gives no contribution to this tail, which 
comes accordingly only from the contribution of T(K, CI). 

For this contribution the dependence of E(k, w) on frequency comes explicitely from the Go term, and it has poles 
for uj = CI — eK-k + Here f2 runs over the frequencies of the singularities of T(K, but in practice the Bose 
factor limits them to some finite range fl < T. This implies in particular that K is bounded (just as above the 
terms produced by the tail of the Bose factor will necessarily have an exponential factor and do not contribute to the 
algebraic tail of n/c). Hence, for large k — > oo, wc obtain from the Go term a single pole located atw« — k 2 /2m. The 
corresponding residue is, from Eq. (|Bip : 



1 f fl 1 f°° Q 

R ^ = -^-H / dnb(-)T(K,fl) = — £ / dn&(-)Imr(K J fi + «0 (B3) 

where Cr is a clockwise contour enclosing only the singularities of r(K, CI) and the last expression (with e — > 0+) is 
the real axis integral obtained by calculating the contour integral from the jump of the imaginary part of T across 
the real frequency axis. 

Let us consider now the general situation where no approximation is made. Actually, provided we replace the bare 
propagator Go by the full propagator G, Eq. (|Bl[) remains valid. We have: 

S(k, iuj„) = TJ2 T(K, k, iuj n )G(K - k, iu> v - iu> n ) = ^- f dQ fc(^) ^ 7 (K, fi; k, iuj„)G(K - k, Q - i ^„)(B4) 

K,i/ l7T ^ C K 

Indeed, considering the self-energy of, say, the f particle, there will be in any diagram a first interaction with a 
| particle, hence correspondingly a G| factor. Merely isolating this factor and gathering the rest into 7 produces 
Eq. (|B4j) . We may express 7 in terms of the standard [l(| full vertex f by isolating the Hartree term in the self-energy. 
This gives, for a contact interaction: 

j(K;k) = g - gT^G^k^G^K - k')t{k' ,K - k';k,K - k) (B5) 
fe' 

where wc have used a four- vector notation k = (k, iuj n ). Here the coupling constant g is related to the scattering 
length a and to the cut-off k c , related to the contact interaction, by g^ 1 = m r /[2/Ka) — ^ " 2m r /k 2 . 

Physically 7(K, iuj v \ k, iuj n ) describes the scattering of two particles k and K — fc, and accordingly it will have the 
same qualitative properties as we have discussed above for T(K, CI) (where we have been quite general for purpose). 
We consider the contributions in Eq. (|B4[) from the singularities of G. Either K is small compared to k, so K — k is 
large and we can replace G by Go, and as above wc will have an exponential factor coming from the Bose factor. Or 
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K is comparable to k, in which case the corresponding singularities of 7 will start around K 2 /Am and the resulting 
singularities for £(k, uj) will be at least u> « k 2 /Am, which makes them irrelevant for the algebraic tail of n^. 

Hence we have only to consider the contributions of the singularities of 7 in the calculation of Eq. (|B4|) . Again those 
with large K (i.e. with small K — k) will be located at high frequency fi, and the Bose factor produces an exponential 
factor which makes them irrelevant. We are left only with contributions arising for K small compared to k, in which 
case G can again be replaced by Go- This produces for S(k, oS) an explicit pole at uj w —k 2 /2m, with a residue: 



Actually it seems physically reasonable that, in the large k limit. 7(K, il; k, ~k 2 /2m) depends only on (K, SI). Indeed 
this limit corresponds to a short range scattering between two | and | particles. In this case, in a way analogous 
to the one discussed in section |n] for the energy formula, this scattering should be essentially described by two-body 
physics. The situation is then similar to the one found in the ladder approximation. The k dependence is just the 
one which is already explicit, and we are left with a constant for the residue, given by a formula like Eq. (|B3[) with a 
modified r(K,0). 
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